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An examination of proton charge radius extractions from e-p scattering data 


John Arrington 

Physics Division, Argonne National Laboratory, Argonne, IL 60439 

A detailed examination of issues associated with proton radius extractions from elastic 
electron-proton scattering experiments is presented. Sources of systematic uncertainty and 
model dependence in the extractions are discussed, with an emphasis on how these may 
impact the proton charge and magnetic radii. A comparison of recent Mainz data to previous 
world data is presented, highlighting the difference in treatment of systematic uncertainties 
as well as tension between different data sets. We find several issues that suggest that larger 
uncertainties than previously quoted may be appropriate, but do not find any corrections 
which would resolve fhe proton radius puzzle. 

PACS numbers; 13.40.Gp,13.40.Gp,14.20.Dh,25.30.Bf 


I. INTRODUCTION 


Five years after the initial extraetion of the proton radius from muonie hydrogen yj], the “proton 
radius puzzle” persists. Measurements based on muonie hydrogen transitions yfl and those based 
on eleetron transitions |3] or eleetron seattering measurements |4-^ disagree at the la level, with 
muonic hydrogen results yielding a radius near 0.84 fm and electron-based measurements yielding 
te ~ 0.88fm, as summarized in Fig.[TJ In light of this, a careful examination of the details of these 
extractions is clearly warranted. Here, we discuss several issues relevant to determining the proton 
radius from electron scattering data. 

In examinin g e xtractions from electron scattering data, we examine the Mainz data and 
global analyses 1^ [SD of world data (excluding Mainz) separately. This is done because the Mainz 
data presents the uncertainties in the data in a significantly different way from most other ex¬ 
periments, making it difficult to perform a meaningful combined analysis. It is also beneficial 
to perform independent analyses to examine consistency between the Mainz data and other mea¬ 
surements at the cross sections level, which can be overlooked in a combined analysis. We also 
discuss some preliminary results from a detailed examination of both Mainz and world data [SI]. 
There are several issues that suggest that larger uncertainties than quoted in previous works are 



















2 



6.7 o 



0 A A 0 0 


«- 

-^ 



pp2013 




electron avg. 



-•- 


scatt. JLab 

pp 2010 



— • — 

scatt. Mainz 

I I I I I I I I 

I I I I I I I I I I I I I I I 

-« 

— 

^ H spectroscopy 


0.83 0.84 0.85 0.86 0.87 0.88 0.89 0.9 


Proton charge radius R [fm] 

ch 


FIG. 1: Extractions of the proton charge radius from muonic iwdrogen measurements 111,121!, hydrogen 
spectroscopy |j^, electron scattering measurements at Mainz 10, Q], and a global analysis of earlier world 
data 10]. The direct average shown is compared to the CODATA-2010 evaluation |0]. Figure courtesy of 
Randolf Pohl. 


warranted. While none of these appear likely to resolve to the discrepancy with muonic hydrogen 
measurements, some issues remain which deserve more detailed examination. 


II. GENERAL ISSUES IN THE EXTRACTION OF THE RADII 


One obtains the charge and magnetic form factors, and from unpolarized 

cross section measurements by performing a Rosenbluth separation IgD which uses the angle- 
dependence at fixed to separate the charge and magnetic contributions. The cross section at 
fixed is proportional to the ’reduced’ cross section aR = tG\j + eG%, where r = Q'^/ (4Mp) 
and = [1 + 2(1 + r) tan^(6*/2)]. At low Q^, the magnetic contribution is strongly suppressed 

except for very small e values, corresponding to large scattering angle. Because of the difficul¬ 
ties in making very large angle scattering measurements at low Q^, a significant extrapolation to 
£ = 0 is required and even sub-percent uncertainties on the cross sections can yield significant 
uncertainties on small contribution from Gm(Q^). 

Because one often combines data from many experiments, each of which has an uncertainty in 
its normalization uncertainty, the normalizations factors of the limited number of large-angle data 
sets have a great impact on the extraction of Gm- If these normalization factors are allowed to 
vary in the fit, which is the most common approach, then a small shift in normalization between 
large and small angle data sets can yield a significant shift of strength between Gp, and Gm over 
a range in values. Polarization observables are sensitive to the ratio Gr/Gm [llOiillJ] and can 
thus provide not only direct information on the form factors, but also improve the determination 
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of the relative normalization of different measurements. Of partieular interest are data sets at low 


values y, 


12 


1411 whieh provide improved extraetions of Gm and additional eonstraints on the 


experimental normalizations. 

Extraetion of the eharge radius from eleetron seattering requires parameterizing the eross sec¬ 
tions to obtain the slope of the form factor at = 0. Many naive extractions use fit functions 
which do not provide sufficient flexibility to accurately describe the low-Q^ data, and often do not 
attempt to estimate the uncertainty associated with the choice of functional form used to fit the 
data. For example, several early extractions were based on linear fits to low data. Such fits will 
always give an underestimate of the radius, based on the observed positive curvature of the form 
factors at low 0^. One would have to have extremely precise data at very low for a linear fit to 
be sufficient 115, 1^ . 

Simply increasing the number of parameters in a simple Taylor expansion or similar fit function 
provides greater flexibility. However, it also leads to increased uncertainty in the extracted radius 
as correlations between different parameters in the expansion allow the impact of variations in one 
term of the expansion to be balanced by changes in other terms. This yields a rapid increase in the 
extracted radius uncertainty as the number of parameters is increased. This may lead to a situation 
where there is no region in which there are enough parameters to accurately reproduce the data 
while still yielding an uncertainty small enough to provide a useful radius extraction [SI]. Such 
analyses must find a balance between fit flexibility and radius uncertainty and, ideally, attempt 
to estimate the error made when truncating the fit function. All of the extractions that we will 

nn 

review in detail ||41-[^ examine the model dependence associated with the functional form used to 
parameterize the data and include at least some estimate of the associated uncertainty. 

The factor r suppresses the magnetic contribution as 0, causing the uncertainties on 

Gm{Q‘^) to increase rapidly as seen in Fig.[2l Because the radius extraction is sensitive to the low- 

behavior of Gm and the most precise data are at higher values, it is particularly difficult to 
reliably extract the magnetic radius. In the analyses of world data ll4|,|5l], the fits exclude high-Q^ 
data t^revent these data from influencing the extraction of the slope. For the analysis of the Mainz 
data [^, the data set extends to ~ 1 GeV^, but the bulk of the data below = 0.5 GeV^ and 
a more flexible fit function is used to provide greater flexibility to fit cross section measurements 
at both low and high Q‘^. 

Radiative corrections are another area requiring special attention. The largest contributions 
to the radiative corrections can be calculated in a model-independent way, although there are 
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FIG. 2: Extractions of Gm and their uncertainties from direct Rosenbluth separations for the Mainz data 
(crosses) and from a global (pre-Mainz) analysis jlTI] (circles). 


exchange (TPE) contributions 1122 . 


small variations between different prescriptions 1118142ih . Other terms, in particular the two-photon 


230, are model dependent as it is necessary to account for the 


possible hadronic states in between the two exchanged photons. The world data analyses y, l50 


include two-photon exchange corrections based on a calculation in a 


ladronic basis including 


only intermediate proton state ll24ll . although estimates of excited states 11251.12611 suggest that their 
contribution is very small at the relevant values. The radiative correction uncertainties quoted 
by the experiments used in these global analyses were typically 1-1.5%, and are assumed to be 


su 


ficient after applying the calculate TPE corrections. We note that the data of Simon, et al. 1127 . 


2811 did not include any uncertainty for radiative corrections and thus tend to have an artificially 


enhanced impact on extractions of the form factors and radius. In the analysis of Ref. [|40, and 
additional radiative correction uncertainty was applied to the Simon data. 

n 

The primary result from the Mainz experiment applies TPE corrections derived for a point 
target [l29n (the “Eeshbach” correction). This correction is exact for = 0 but has no depen¬ 
dence. Because the radius is the slope of the form factors at = 0, it seems unlikely that a 
-independent correction will be sufficient. The model-dependent TPE calculations mentioned 
above agree with the Eeshbach correction at = 0, but as increases they tend to decrease, 
going to zero before changing sign and growing in magnitude above ~ 0.3 GeV^. There are 
several TPE calculations meant to be appropriate at low [fl Q IsolH. and they are all in 
good agreement at low as shown in Ref. [l35ll . Very recently, this change of sign relative to 
the = 0 limit was confirmed by comparisons of electron-proton and positron-proton scattering 
for Q2 ~ 1 and 1.5 GeV^ H 1^ . This supports the idea that the = 0 calculation is not 
appropriate and a more complete TPE correction is required. 
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The question of TPE eorreetions in the Mainz data was first examined in Refs. OSL 


M- 


Ref. 1I39I1 shows a direet eomparison of the extraeted value of UpGe/Gm with and without TPE 
eorreetions from Ref. |32], whieh are expeeted to be valid up to ~ 0.1 GeV^. As noted in Q, 
the eorreetion on UpGe/Gm is relatively small, below 1%. However, this eorreetion is larger than 
the linear sum of the statistieal, systematie, and model uneertainties. It is, therefore, a eritieal 
eorreetion for an extraetion aimed at sueh hig^reeision, and elearly neeessary for a preeise ex- 
traetion of the eharge and magnetie radii. Ref. does not inelude any uneertainty assoeiated with 
TPE eorreetions, but does inelude an extraetion of the radius after applying hadronie eorreetions 


with the proton intermediate state ll24ll . The ehange in the eharge radius is 0.004 fm, roughly 1/3 


of the total uneertainty [TI], while the magnetie radius ehanges by 0.022 fm, more than the total 
quoted uneertainty. 


III. EXAMINATION OF THE MAINZ ANALYSIS 


As noted earlier, the extraetion of the uneertainties as well as the breakdown of different types 
of uneertainties in the reeent Mainz data set is signifieantly different from other experiments. We 
deseribe the approaeh used in Ref. and then diseuss potential implieations on the uneertainties 
of the extraeted radii in the Mainz analysis, as well as independent fits to the Mainz eross seetion 
data. 


A, Uncorrelated systematic uncertainty 

The uneorrelated systematie uneertainties were determined by performing a fit to the full data 
set using only the pure eounting statisties for uneertainties. The differenee between the data and fit 
for eaeh subset (eaeh independent energy-speetrometer eombination) was examined, and a sealing 
faetor was determined for eaeh data set whieh, when applied as a seale faetor enhaneement to the 
uneertainties from the eounting statisties on every data point, yielded a seatter that was approx¬ 
imately eonsistent with the enhaneed statistieal uneertainty. The goal is to provide a redueed 
value eloser to unity, with xl ~ 1-14 for the final Mainz fit to the eross seetions with the sealing 
faetors applied. This proeedure yields the minimal uneertainty neeessary to aeeount for the non- 
statistieal seatter of the data, but is insensitive to any sourees of error whieh may be eorrelated 
with the kinematies or operating eonditions of the experiment, e.g. beam energy or speetrometer 
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angle offsets, approximations in the radiative correction procedures, or subtraction of target cell 
wall contributions. In fact, because the final reduced chi-squared is still above one, the final un¬ 
correlated systematic uncertainty is somewhat below the minimum necessary to account for the 
observed scatter. 

This rescaling procedure is relatively unusual; nearly all other experiment made direct esti¬ 
mates of uncertainties or upper limits for various sources of uncertainty which may be treated as 
uncorrelated in the fit. This uncorrelated systematic is determined and added in quadrature to the 
statistical uncertainty. If we convert the Mainz scaling factors to independent systematic uncer¬ 
tainties using this standard approach, they correspond to uncertainties that average 0.25%, but vary 
from 0.02% to 2% with the smallest systematic uncertainties generally being applied to the data 
with smallest statistical uncertainties. 

B, Correlated systematic uncertainty 

Most experiments provide relatively small data sets, typically tens of cross section measure¬ 
ments covering a range of and e values. For such data sets, correlated errors, e.g. associated 
with kinematic-dependent corrections, can be well represented by applying an additional uncor¬ 
related uncertainty to each point. A modest 0.5% contribution to the uncertainty on each cross 
section provides flexibility to cover an arbitrary correlated uncertainty at the few tenths of a per¬ 
cent level. For the Mainz measurement there are 1422 cross sections, 10-50 times more than 
most experiments, so any limited kinematic region will have many more data points, reducing 
the impact of the uncorrelated uncertainty by the square root of the number of points. Thus, try¬ 
ing to represent small correlated effects as uncorrelated uncertainties would require much larger 
contributions. 

This effect is made worse by the fact that of the 1422 data points, there are only 638 indepen¬ 
dent kinematic settings. In several cases, multiple repeated measurements were taken at the same 
kinematic setting, one after the other. For the given procedure - inflating the counting statistics by 
a scaling factor intended to yield a reasonable chi-squared for each data set - it doesn’t matter that 
there are multiple repeated measurements in the data set. However, using the more conventional 
approach of applying a fixed systematic uncertainty to each point, a set of N repeated measure¬ 
ments would artificially reduce the impact of the systematic uncertainty by a factor of \/]V for 
this kinematic setting. If the data are rebinned into their 638 independent points and a systematic 
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FIG. 3: Size of the correlated systematic parameter a for Spectrometer A (diamond), B (triangle), and C 
(circle) for each beam energy setting of Ref. 


applied that yields a reduced chi-squared value near unity, the uncertainties tend to increase more 
where the scaled uncertainties were very small, i.e. the very high statistics points or the kinematics 
with small scaling factors. This ends up increasing the low data uncertainties more and yields 
a larger uncertainty on the extracted radius [SD. 

Because of the limitations of including only uncorrelated systematic uncertainties on such a 
large data set, the A1 collaboration treated correlated systematic uncertainties independently. They 
separated the full data set into 18 subsets, each corresponding to a single spectrometer and fixed 
beam energy. They then simultaneously applied a correction factor, proportional to the scattering 
angle, to all points with each of these 18 subgroups and refit the data. The correction varied 
from 0% at the smallest angle to a setting-dependent maximum value at the largest angle for each 
setting. The size of this maximum correction, the parameter “a” from Ref. |^, is shown for all 18 
spectrometer/energy combinations in Fig.[3l with most settings having a value between 0.1% and 
0.25%. In addition, that there is also a separate correlated systematic which accounts for variation 
of the cross section with the elastic tail cut, which is evaluated separately and then combined with 
correlated systematic mentioned above. As this is the smaller contribution, we focus here on the 
correction that is taken to be linear in 9. 

Note the in the supplemental material of Ref. [^, the systematic correction does not go to zero 
for the smallest scattering angle of each subset, but goes to zero at 6^ = 0. If the normalization 
factors of the data subgroups are allowed to vary, the two procedures are equivalent. However, if 
the normalizations are not allowed to vary, applying the corrections factors as published will yield 
a much larger correlated error. 
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C. Normalization uncertainty 

The combination of the uncorrelated and correlated uncertainties applied to the Mainz data are 
extremely small compared to other measurements and, by themselves, represent an incomplete 
estimate of the experimental uncertainty. The data are broken up into 34 subgroups, each of which 
has an independent, unconstrained normalization factor. Thus, the uncorrelated and correlated 
systematic uncertainties described above need only account for the variation of any corrections 
over the kinematics of the individual subgroups which consist of between 18 and 68 independent 
kinematics each (treating multiple runs at identical kinematics as single points), as the overall 
normalization factor will account for any average correction. Note, however, that there are 34 
different normalization subgroups while the correlated systematic uncertainties are applied over 
the 18 independent beam energy-spectrometer combinations. It would seem more consistent for 
apply the correlated systematic over each of the 34 normalization subgroups, and the potential 
impact of this, as well as the choice of functional form for the correlated systematic, will be 
discussed in the following section. 


D. Treatment of the uncertainties 


Because the uncertainties of the measurement are separated into uncorrelated, correlated, and 
normalization uncertainties, a proper evaluation of the radius uncertainties must account for all 
of these. Fits which take the Mainz cross section data as quoted and do not allow the normal¬ 
ization factors to vary, e.g. as in Ref. [l40n . will yield artificially small uncertainties in the radius 
extractions, as discussed in Ref. [SD. 

In addition, questions have been raised about missing contributions to the uncertainties. No 
uncertainties associated with TPE corrections are included, and there are additional uncertain¬ 
ties associated with approximations made in radiative correction procedures. These are neglected 
in dbl], based on the assumption that these will be contained in the small correlated systematic 
uncertainties applied for other effects, but no argument is made to support these uncertainties be¬ 
ing negligible compared to the typical 0.1-0.2% correlated systematic (Fig. |3]). They also quote 
uncertainties on the knowledge of the beam energy and spectrometer scattering angles, but do not 
account for these in the cross section uncertainties. While the impact of these kinematic uncer¬ 
tainty on the cross sections is generally very small, the corrections for both energy and angle can 
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be as large as 0.2%, mainly at low energy and small scattering angle, and are strongly kinematic 
dependent. Thus, it is not clear that they should be neglected in comparison to statistical uncer¬ 
tainties at the 0.2% level and correlated systematics which are as low as 0.1% for some settings. 

It has also been noted ISl] that the target cell wall contributions, subtracted from the data based 
on a calculated spectrum, don’t match the observed contribution. Based on the visible difference 
between the data and simulated spectrum in the region between the nuclear elastic and e-p elastic 
peaks (Fig. 8b of Ref. [a), it was estimated that this subtraction underestimates the cell con¬ 
tributions by 1.2% [qI]. This is a large effect, both compared to the quoted uncertainties of the 


measurement and to the size of the subtraction (below 4% for most settings OQIl l. While this is the 


largest contribution to the correlated systematic uncertainty, it is still taken to be roughly 0.1-0.2% 
for nearly all settings. Any overall normalization error caused by an underestimate of the back¬ 
ground subtraction will be removed in the fit, but with the subtraction varying from a few percent 
to nearly 10% for spectrometer B at forward angles ll39ll . it’s not clear that the kinematic variation 
over each of the data subsets can be constrained to be below the 0.1-0.2% level. 

In addition to noting that the correlated systematic uncertainties applied are very small and may 
neglect important contributions, we also note that the impact of these uncertainties on the extracted 
form factors and radii are evaluated within a single model; a shift of the cross sections which is 
linear in 6 over each of the 18 spectrometer-beam energy combinations. Taking the correction to 
scale with quantities other than 6 can yield larger or smaller corrections, with an increase 50% or 
more in some models where the correction is not linear in angle [SI] . Note that the procedure always 
yields a fixed correction between the smallest and largest angle settings of a given data subset, so 
it is only the form of the variation over this subset that is changed in these tests. The impact of 
the correlated systematic uncertainty is also significantly increased if it is taken as a correction 
to each data set with independent normalization, which seems a more consistent approach given 
the breakdown of uncertainties into uncorrelated, correlated, and normalization. Given the issues 
with the size of the correlated systematics noted above and the model dependence of converting 
these systematic effects into an uncertainty on the extracted radius, a conservative approach would 
appear to yield significantly larger radius uncertainties associated with the correlated systematics. 

Note that the detailed comparisons in Ref. [|8l] are performed for the rebinned version of the 
Mainz data. When examining the original analysis procedure [^, the correlated systematic uncer¬ 
tainty on the extracted radii is more sensitive to both the functional form chosen to represent the 
correlated corrections and to the question of whether the correction is applied to each normaliza- 
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tion subgroup or each beam-spectrometer combination. 


IV. FITTING WITH BOUNDED z EXPANSION 


As noted earlier, it is important that the fit function be flexible enough to adequately reproduce 
the data without being so flexible that it does not provide a meaningful constraint on the extracted 
radii. This is especially critical for the magnetic radius, where the precise high-Q^ data can in¬ 
fluence the fit more than the low-Q^ data which are directly sensitive to the radius. Typically, a 
fit function is selected and the number of parameters is chosen to be as small as possible while 
still providing a reduced chi-squared that is close to the minimum value obtained for many pa¬ 
rameters. In some cases, fit functions designed to help minimize the impact of the high-Q^ data 
on the low-Q^ fit are chosen, e.g. spline functions [Q] or continued fraction fits [15]. However, it 
is still difficult to determine how many parameters are sufficient for a reliable fit. Figure 9.21 of 
Ref. [Q] shows the extracted charge and magnetic radii vs. the number of parameters for a variety 
of different functional forms. For the charge radius, the extracted radii are relatively consistent for 
fits with 10-15 parameters, although there is a significant spread (~ 0.02 fm) between the values 
using different fit functions. For more than 16 parameters, the spline fits yield much smaller radii, 
but this is presumably in the region where the uncertainties become very large and so the shift of 
the central value may not be outside of the fit uncertainty. For the magnetic radius, the situation 
is noticeably worse. There is a narrow window, 6-8 parameters, where the radii are relatively 
stable and then by 10 parameters the different functional forms yield radius values differing by 
nearly 0.2 fm. As noted in Ref. [|42|] . some of this erratic behavior is associated with fits that have 
unphysical behavior, e.g. poles in the form factor and oscillatory behavior of the proton charge 
density, p(r) at very large values of r. 


One approach to this problem was presented in Refs, oi |42l], where a parameterization was 


chosen that constrains the large-r behavior of the form factors. By including physical constraints 
on the behavior of the form factors, one can avoid the possibility that insufficiently flexible fit 
functions may yield a poor fit radius to better reproduce high-Q^ data. 

Another approach that can help address issues of over- or under-fitting data is the use of the 
bounded 2 ; expansion [l43n . In this method, the form factors are parameterized as a polynomial in 
2 ; rather than where 


tcuty to) 


\/^CUt t V^fcut ^ 

\/fcut t -f \/fcut to 


( 1 ) 
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where tcut = 4m^ and is a free parameter. The true form faetor is guaranteed to be in the space 
of this polynomial expansion for large enough number of parameters, and sum rules exist which 


limit the size of the coefficients in the expansion. The details o 


the proton charge and magnetic form factors in Refs. [S, 


44|-l46|. 


the 2 ; expansion are discussed for 


The constraints on the coefficients allow us to estimate bounds which can be applied to the 
individual fit coefficients. Applying such a bound to the fit prevents the uncertainties from growing 
out of control as more parameters are added, because it damps the oscillatory behavior that can 
occur if each term cancels the error from the previous term. The bounded 2 : expansion provides 
fits which are more stable in both the extracted radius and uncertainty as one increases the number 


of fit parameters [Si |44|, 1451]. This allows the number of parameters to be large enough that the 


fit is not limited by the truncation of the Taylor series, while still providing a meaningful, though 
perhaps larger, uncertainty on the radii. It also helps to decouple the parameters needed to fit 
the low data from those important at high making the fits somewhat more robust against 
potential error in measurements at larger values where the form factor has little sensitivity to 


44 


45]], and if the bound 


the radius. The estimate of the bounds, however, is model-dependent [Si 
applied it too tight, it can bias the extraction of the radius. 

The ability to bound the fit coefficients and go to high order fits is particularly useful when 
comparing the analysis of Mainz data and the world data set, as it allows for the same functional 
form and number of parameters to be used in both cases. The charge radii extracted from the 


Mainz data 0] and world data [141,15l] are in good agreement, but the magnetic radii are signifi¬ 
cantly different. The fact that these analyses use very different fit functions and have to select a 
range of parameters tailored to the size and precision of the data sets makes it difficult to deter¬ 
mine the role of the model dependence of the fits. With the 2 : expansion, different data sets, as 
well as different ranges of a single data set, can be examined in a way that minimizes model 
dependence associated with choosing different fit functions. A detailed analysis discussing the re¬ 
maining model dependence associated with the 2 : expansion fits and comparing consistent analyses 
of Mainz and world data has been undertaken [SI]. This comparison shows that the charge radius 
from both Mainz and world data are still inconsistent with muonic hydrogen results, and that the 
discrepancy in magnetic radii extracted from the Mainz and world data persists. 

Other analyses have used functions with fewer parameters, requiring only that the number is 
sufficient to provide an approximate plateau in the chi-squared value of the fit and the extracted 
radii. This is done because the uncertainties grow significantly with the number of parameters. 
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and so fits are generally taken with the minimum number of parameters required to reasonably 
fit the data. For the 2 ; expansion, the results are independent of the number of parameters onee 
this beeomes large, and so the most reliable approaeh is to have several more parameters than 
is neeessary to obtain a reasonable ehi-squared value, to avoid under-fitting of the data. The 
fit uneertainties thus tend to be somewhat larger than quoted in previous results, espeeially when 
taking eonservative estimates of the eoeffieient bounds. However, the fits should yield more robust 
estimates of the uneertainties as they avoid the large dependenee of the uneertainty on the number 
of parameters used in the fit. 


V. CONSISTENCY BETWEEN MAINZ AND WORLD DATA 


As noted earlierj_^there are signifieant tensions between the Mainz data set and world data. 
Figure 10 of Ref. eompares the Mainz fit to previous world data, and shows a signifieant 
disagreement between the Mainz measurement and nearly all low-Q^ extraetions of Gm- The 
disagreement in Gm is roughly 3-4%, eorresponding to an 6-8% eross seetion differenee if this 
were explained entirely by normalization faetors, as suggested in Bfil]. However, while the Mainz 
analysis yields values of Gm that are several pereent above world data, the Mainz results for Ge 
are systematieally below world data for 0.2 < < 0.8 GeV^. Thus, a simple normalization 

eorreetion eannot resolve the diserepaney. Note that the world data results do not have TPE eor- 
reetions applied, but at these values the eorreetions are relatively small and for a eomparison 
to the Mainz result, the Feshbaeh eorreetion should be used. Applying this eorreetion to the world 
data would have a small effeet that would deerease Gm and inerease Ge, further inereasing the 
tension with the Mainz analysis. A reeent analysis [l47ll examined both Mainz and world eross 
seetion data, and extraeted the TPE eontribution using a phenomenologieal approaeh. While the 
exaet TPE extraeted at these low values may not be determined preeisely, this proeedure ap¬ 
plies a eonsistent eorreetion to both Mainz and world data and a signifieant tension is observed, in 
partieular in the magnetie form faetor. 

One ean also see the disagreement in Eigure 19 of dg], whieh shows the normalization faetors 


applied to previous measurements as determined from a global fit. All of the world data sets shown 
require an inerease in their normalization, with roughly half of these renormalized by 4% or more. 
This ineludes several data sets whieh are shifted by 2-3 times their quoted normalization uneer- 
tainties. As noted earlier, the differenee in the way the uneertainties are separated in the Mainz 
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data will lead to it receiving an artificially high weight in the fit, but this analysis provides a di¬ 
rect measure of the large relative renormalization factors required to improve consistency between 
Mainz and world data. 


VI. CONCLUSIONS 


The large number of high-precision data points in the recent Mainz experiment by the A1 
collaboration requires special treatment of the uncertainties. The release of the data includes un¬ 
correlated uncertainties (statistical and systematic), correlated systematic uncertainties, and nor¬ 
malization factors for 34 different subsets of the data. This leads to two concerns with analyzing 
these data. First, it is extremely important to account for all of these uncertainties, in particular the 


normalization factors and corre 
these data, as done in [y, [sl. 


ated systematic uncertainties, in any extraction of the radius from 


4611 . Second, a simple global analysis of the Mainz data with other 


cross section or polarization observable measurements will give too much weight to the Mainz 
data, as the uncertainties given for each cross section point represent only a small fraction of the 
total uncertainty. 

Given recent measurements supporting the importance of the hadron structure-dependent TPE 


corrections Obi. 


3711 . it is clear that a -independent TPE correction is not sufficient and that an 


uncertainty associated with the model dependence of the TPE correction needs to be included. We 
note this and other contributions that are not included in the evaluation of the correlated systematic 
uncertainties, but which can have a significant impact on the uncertainty of the extracted radius. 
In addition, different assumptions about the kinematic dependence of these unknown systematic 
corrections can noticeably increase the impact of these corrections on the extracted radii. Between 
missing contributions to the total systematic uncertainties and the model dependence of evaluating 
the impact of these corrections, it appears that the systematic uncertainties associated with the 
radius extraction of are likely to be significant underestimates of the true uncertainty. 

Evaluating the model dependence of such fits is important, and while one can select fit functions 
and ranges of parameters which appear to yield good fits to the data with reasonable uncertainties, 
it is difficult to cleanly determine if one is under-fitting or overfitting the data, either of which can 
significantly modify the extracted radius and uncertainties. We argue for the use of the bounded 
2 ; expansion, which allows a large number of parameters without concern about overfitting and 
without the dramatic loss of precision that comes with unbounded fits at high order. This procedure 
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yields somewhat larger uncertainties, but yields a more reliable determination of the uncertainty 
in the extracted radius. 

Initial results from bounded 2 ; expansion fits to both the Mainz data and world data [[80 yield 
consistent charge radii which are still significantly higher than the muonic hydrogen results. How¬ 
ever, there are clearly inconsistencies between the Mainz data and other world data, both at the 
cross section level and in the extraction of the magnetic radius, where the Mainz data set yields 
a much smaller magnetic radius. Without a better understanding of the origin of the tension be¬ 
tween the different data sets, it is difficult to make a clear and rigorous statement about the present 
uncertainty on the proton’s charge radius as derived from elastic electron scattering. Based on the 
considerations presented here and other examinations of the model dependence of the radius ex¬ 


tractions p. 


421], a recommendation for a radius and uncertainty based on published fits to electron 


scattering data [|4l-l^ is presented in Ref. ll48h . 

This work was supported by the U.S. Department of Energy, Office of Science, Office of Nu¬ 
clear Physics, under contract DE-AC-06CH11357 
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